Importance of bath dynamics for decoherence in spin systems 
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We study the decoherence of two coupled spins that interact with a chaotic spin-bath environment. 
It is shown that connectivity of spins in the bath is of crucial importance for the decoherence of the 
central system. The previously found phenomenon of two-step decoherence (Phys. Rev. Lett. 90, 
210401 (2003)) turns out to be typical for the bath with a slow enough dynamics or no dynamics. 
For a generic random system with chaotic dynamics a conventional exponential relaxation to the 
pointer states takes place. Our results confirm a conjecture of Paz and Zurek (Phys. Rev. Lett. 
82, 5181 (1999)) that for weak enough interactions the pointer states are eigenstates of the central 
system. 
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It is commonly accepted that decoherence by nuclear 
spins is the main obstacle for realization of quantum com- 
putations in magnetic systems; see, e.g., discussions of 
specific silicon [lj and carbon Q based quantum comput- 
ers. Therefore, understanding the decoherence in quan- 
tum spin systems is a subject of numerous works (for 
a review, see Refs. 0, 0). The issue seems to be very 
complicated and despite many efforts, even some basic 
questions about character of the decoherence process are 
unsolved yet. Most of the problems cannot be solved 
analytically, in particular if there is more than one spin 
in the central system, but we can use the computer to 
simulate the dynamics and find useful information. 

An unusual two-step decoherence was reported in 
Ref. |. This is an important phenomenon since it im- 
plies that, generally speaking, the observation of the Rabi 
oscillations does not guarantee access to sectors of the 
Hilbert space that may be essential for efficient quantum 
computation. Its origin is still poorly understood; it was 
described analytically in a framework of an exactly solv- 
able model of noninteracting spins in the bath Q but it 
is not clear how sensitive it is to the details of spin-spin 
interactions. In the real world, the environment has its 
own dynamics, which could be much slower or compara- 
ble to the central dynamics. First attempts to investigate 
numerically the effects of the environment dynamics Q 
did not lead to definite conclusions. 

The behavior of an open quantum system crucially de- 
pends on the ratio of typical energy differences of the 
central system SE C and the energy E ce which character- 
izes the interaction of the central system with the envi- 
ronment. The case SE C -C E ce has been studied exten- 
sively in relation to the "Schrodinger cat" problem and 
the physics is quite clear d, 0] : As a result of time evo- 
lution, the central system passes to one of the "pointer 
states" 0] which, in this case, are the eigenstates of the 
interaction Hamiltonian. The opposite case, 5E C 3> E ce 
is less well understood. There is a conjecture that in this 



case the pointer states should be eigenstates of the Hamil- 
tonian of the central system, but this is proven only for 
a very simple model [10( . On the other hand, this case is 
of primary interest if, say, the central system consists of 
electron spins whereas the environment are nuclear spins 
(e.g., if one considers the possibility of quantum compu- 



11, 12]). 



tation using molecular magnets 

In fact, as we will show, the selection of an eigenstate 
as the pointer state is also determined by the state and 
the dynamics of the environment. Elsewhere 13. Il4j. we 
have already shown that if the environment is a spin glass 
and initially in the ground state, then independent of 
the initial state of the central system, the central system 
relaxes to a state that is very close to its ground state: 
The ground state is selected as the point state of the 
central system. In this Letter, we consider a realistic 
model of decoherence of a system of two spins by an 
environment of nuclear spins at elevated temperatures. 
We will demonstrate that the decoherence of the central 
system depends in a significant, nonintuitivc manner on 
the details of the dynamics of the environment. 

We consider a generic quantum spin model described 
by the Hamiltonian H = H c + H ce + H e where 
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are the Hamiltonians of the central system, the environ- 
ment, and the interaction between the central system and 
the environment, respectively. In Eq. ([I]), the exchange 
integrals J and Clf*) determine the strength of the in- 
teraction between spins S„ = (S%, S%, S%) of the central 
system H c , and the spins I„ = (I*, in the envi- 
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ronment H e , respectively. The exchange integrals A[°j 
control the interaction H ce of the central system with its 
environment. In Eq. {T]), the sum over a runs over the x, 
y and z components of spin- 1/2 operators S and I. 

In the sequel, we will use the term "Heisenberg-like" 
H e to indicate that each off) is a uniform random num- 
ber in the range [— fi| J\, fi|J|], being a free parameter. 
We will consider two different kinds of H ce , namely ro- 
tational invariant Heisenberg interactions A^ = A and 

(x) (v) 

"Ising-like" interactions for which A 4 - ■ = A\ ■ = and 

A\ J are dichotomic random variables, taking the values 
± A. Obviously, if H ce is "Ising like" , the total magneti- 
zation of the central system (M — Sf+S%) is a conserved 
quantity. 

As we demonstrate in this Letter, the connectivity of 
the spins in the environment affects the decoherence in 
a nontrivial manner. We characterize this connectivity 
by the number K, the number of environment spins with 
which a spin in the environment interacts. If K = 0, 
each spin in the environment interacts with the central 
system only. If K — 2, the structure of the environ- 
ment is assumed to be that of a ring, that is each spin 
in the environment interacts with two other spins only. 
Likewise, K — 4 and K — 6 correspond environments 
in which the spins are placed on a square or triangular 
lattice, respectively and interact with nearest-neighbors 
only. If K — N, each spin in the environment interacts 
with all the other spins in the environment and, to give 
this case a name, we will refer to this case as "spin glass" . 

If the Hamiltonian of the central system H c is a per- 
turbation, relative to the interaction Hamiltonian H ce , 
the pointer states are eigenstates of H ce 0. In the 
opposite case, that is the regime |A| <C \J\ that we 
explore in this paper, th e p ointer states are supposed 
to be eigenstates of H c [lfj. The latter are given by 
|1} = |Tx) = I TT>, |2) = \S) = (| n> - 1 IWV2, 
|3) = |T ) = (| U) + I it»/V2, and |4) = |T_x) = | jj), 
satisfying H C \S) = (3J/4)|5) and H B \fy = (-J/4)|Ti) 
for i = — 1,0, 1. To check this conjecture is one of the 
main aims of our simulations. 

The simulation procedure is as follows. We generate 
a random superposition \(f>) of all the basis states of the 
environment. This state corresponds to the equilibrium 
density matrix of the environment at infinite tempera- 
ture. The spin-up - spin-down state (| f J,)) is taken as 
the initial state of the central system. Thus, the initial 
state of the whole system reads \^f(t = 0))) = |t|) |</>) and 
is a product state of the state of the central system and 
the random state of the environment which, in general is 
a (very complicated) linear combination of the 2 N basis 
states of the environment. In our simulations we take 
N = 16 which, from earlier work 1^, 14|, we know is suf- 
ficiently large for the environment to behave as a "large" 
system. For a given, fixed set of model parameters, the 
time evolution of the whole system is obtained by solving 
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FIG. 1: (Color online) The time evolution of the real part 
of the off-diagonal element P23 (left panel) and the diago- 
nal elements pn, . . . ,^44 (right panel) of the reduced density 
matrix of a central system (with J = —5), interacting with 
a Heisenberg-like environment H e (with fl = 0.15) via an 
isotropic Heisenberg Hamiltonian H ce (with A = —0.075 ) 
for different connectivity numbers K of the spins in the envi- 
ronment: (a) K = 0; (b) K = 2; (c) K = 4; (d) K = 6; (e) 
K — N. The absolute values of the off-diagonal elements p\i , 
pis, P14, P24, and P34 are less than 0.002 (results not shown). 



the time-dependent Schrodinger equation for the many- 
body wave function |\&(t)), describing the central system 
plus the environment. The numerical method that we 
use is described in Ref. [l5|. It conserves the energy of the 
whole system to machine precision. We monitor the ef- 
fects of the decoherence by plotting the time dependence 
of the matrix elements of the reduced density matrix of 
the central system. As explained earlier, in the regime of 
interest |A| <C \J\, the pointer states are the eigenstates 
of the central systems, hence we compute the matrix el- 
ements of the density matrix in the basis of eigenvectors 
of the central system. 

In Fig.[TJ we show the time evolution of the elements of 
the reduced density matrix of the central system for dif- 
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ferent connectivity numbers K with isotropic Heisenberg 
interactions between the central system and the spins in 
the environment. We conclude that: 
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In agreement with earlier work [B|, Q , we find that 
in the absence of interactions between the environ- 
ment spins (K = 0) and after the initial fast decay, 
the central system exhibits long-time oscillations 
(see Fig. QJ a) (left)). In this case and in the limit of 
a large environment, we have Q 



Re p23 (t) 



1 1 - bt 2 )e- ct 
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6 3 



cos tut, 



(2) 



where b = NA 2 /4, c = 6/2 and lu = J - A. Equa- 
tion ^ clearly shows the two-step process, that is, 
after the initial Gaussian decay of the amplitude 
of the oscillations, the oscillations revive and their 
amplitude levels of [(| . Notice that because of con- 
servation laws, this behavior does not change if the 
environment is described by an isotropic Heisen- 
berg Hamiltonian = f2 for all a^J, and j), 



whatever the value of K. If, as in Ref. [5J, we take 
A,H = A^] = A^ z ) £ [0, Al random instead of the 
same, the amplitude of the long-living oscillations 
is no longer constant but decays very slowly [5j . 

2. The presence of Heisenberg-like (non-isotropic, ran- 
dom) interactions between the spins of the environ- 
ment leads to a reduction and a decay of the am- 
plitude of the long-living oscillations (see Fig. QJb- 
e)(left)). The larger the connectivity number K, 
the faster is the decay of the amplitude. When K 
reaches its maximum K = N (spin glass), the sec- 
ond step in the decohcrence process is no longer 
separated from the initial decay. In fact, it seems 
as if it has merged with the final stage of the first 
step (see Fig. []Je) (left)). For K = N, the time evo- 
lution of P23 (t) can be fitted well by the function 



Re y023 (t) 
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(1 - b't 2 )e 



cosw't, (3) 



with a' = 0.13403, b' = 0.00659, c' = 0.01085 
and ui' — 5.01037. For comparison, the values 
that enter Eq. © arc b = 0.0225, c = 0.01125 
and uj — 4.925. It is of interest to note that 
if the dynamics of the environment is sufficiently 
slow (Q w 0.01 in Ref. [Ej]), this dynamics appar- 
ently does not affect the decoherence of the central 
spins Thus, the effectiveness of the decoherence 
process not only depends on K but also on the de- 
tails of the interactions within the environment. 

3. According to the general picture of decoherence 0, 
for an environment with nontrivial internal dynam- 
ics that initially is in a random superposition of all 
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FIG. 2: The time evolution of the real part of the off-diagonal 
element P23 of the reduced density matrix of a central system 
(with J = —1), interacting with a Heisenberg-like environ- 
ment H e (with Q = 0.15) via an Ising-like Hamiltonian H ce 
(with A = 0.075 ) for different connectivity numbers K of the 
spins in the environment: (a) K = 0; (b) K — 2; (c) K = 4; 
(d) K = 6; (e) K = N. 



its eigenstates, we expect that the central system 
will evolve to a stable mixture of its eigenstates. 
In other words, the decoherence will cause all the 
off-diagonal elements of the reduced density matrix 
to vanish with time. Moreover, the weight of the 
degenerate eigenstates \T ), |T_i), and |Xi) in this 
mixed state are expected to be the same. As shown 
in Fig. [TJb-e) (right), our simulations confirm that 
this picture is correct in all respects. Furthermore, 
the results depicted in Fig. [TJb-e) (right) demon- 
strate that the connectivity number K has no effect 
on the value of pn, P22, P33 and P44 for long times. 



Next, we change the interactions between central sys- 
tem and the spins in the environment from Heisenberg- 
like to Ising-like but keep the interactions between dif- 
ferent spins in the environment Heisenberg-like. In fact, 
in our simulations, the Hamiltonian H e is the same for 
both cases. As explained earlier, in the Ising-like case, 
the total magnetization of the central system is conserved 
during the time evolution. Thus, as the initial state of 
the central system is (\S) + |Xb))/v2, at any time t the 
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TABLE I: Frequency of oscillation to" and decoherence time 
t for different connectivity K. In case 1 and case 2, for each 
K, the values of fl-"' are the same but the values of A,-" are 
different. In case 2 and case 3, for each K, the values of 
are the same but the values of flj) are different. 



Case 




K= 2 


K = 4 


K = 6 


K = N 




1 


n 

UJ 

n 

T 


1.015 
212.9 


1.015 
235.3 


1.016 
302.9 


1.017 
447.4 


2 


u" 

ft 

T 


1.017 
93.07 


1.021 
105.0 


1.022 
110.4 


1.027 
182.1 


3 


ft 

UJ 

ft 

T 


1.021 
106.3 


1.022 
110.0 


1.022 
111.7 


1.027 
173.2 



central system and environment is much smaller than 
the coupling between the spins in the central system, the 
pointer states are the eigenstates of the central system. 
Both these observations arc in concert with the general 
picture of decoherence [§]. 

Furthermore, we have demonstrated that, in the case 
that the environment is a spin system, the details of this 
spin system are important for the decoherence of the cen- 
tral system. In particular, we have shown that (1) chang- 
ing the internal dynamics of the environment may change 
the qualitative features of the decoherence of the central 
spin system, and that (2) the dependence of the decoher- 
ence time of the central spin system on the connectivity 
of the interactions between spins of the environment is 
counterintuitive. 



state of the whole system can be written as 

|*(t)) = |2)|^(i)) + |3)|0a(t)), (4) 

where |02 (*)} and |<fe(*)) denote the states of the envi- 
ronment. In other words, only P22(t), P23(t)> p32(t), and 
P33(t) can be nonzero. 

On general grounds, we may expect that the presence 
of an additional conservation law slows down the deco- 
herence and indeed, as shown in Fig. [21 this is the case. 
Note that the results of Fig. [2] have been obtained for 
J = — 1 instead of for J = —5, the value used to compute 
the results shown in Fig. [1] (the latter value was chosen 
to facilitate the comparison with the results of Ref. [H), 
but this factor of five in the value of J cannot account for 
the large difference in the observed decoherence times. 

From Fig. [21 we conclude that 

1. We never observe the two-step process that we find 
in the case of Hciscnberg-like H ce . For K — (see 
Fig. Ufa)), there is no decoherence. 

2. For K > 0, Re P23{t) vanishes with time, in agree- 
ment with the general picture of decoherence 
However, quite unexpectedly, the rate of decoher- 
ence increases with K, in contrast to the case of 
Heisenberg H ce in which the rate of decoherence 
decreases with K. The data presented in Figs.[2fb- 
d) can all be fitted very well by 

Re P23 (t) = V*/ r "cosw"i, (5) 

where u>" » |J| and the values of t" depend on K. 
In Table [H we give some typical results for these pa- 
rameters. Cases 1 and 2 illustrate that (random) 
changes in H ce may affect the value of the deco- 
herence time t" significantly but the general trend, 
the increase of t" with K seems generic. 
In conclusion, we have shown that (1) the pure quan- 
tum state of the central spin system evolves into the clas- 
sical, mixed state, and (2) if the interaction between the 
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